Nonperturbative effects on the ferromagnetic transition in repulsive Fermi gases 
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It is generally believed that a dilute spin-^ Fermi gas with repulsive interactions can undergo a ferromagnetic 
phase transition to a spin-polarized state at a critical gas parameter {k r a) c . Previous theoretical predictions of 
the ferromagnetic phase transition have been based on the perturbation theory, which treats the gas parameter as 
a small number. On the other hand, Belitz, Kirkpatrick, and Vojta (BKV) have argued that the phase transition 
in clean itinerant ferromagnets is generically of first order at low temperatures, due to the correlation effects 
that lead to a nonanalytic term in the free energy. The second-order perturbation theory predicts a first-order 
phase transition at (k f a) c = 1.054, consistent with the BKV argument. However, since the critical gas parameter 
is expected to be of order 0(1), perturbative predictions may be unreliable. In this paper we study the non- 
perturbative effects on the ferromagnetic phase transition by summing the particle-particle ladder diagrams to 
all orders in the gas parameter. We consider a universal repulsive Fermi gas where the effective range effects 
can be neglected, which can be realized in a two-component Fermi gas of 6 L1 atoms by using a nonadiabatic 
field switch to the upper branch of a Feshbach resonance with a positive s-wave scattering length. Our theory 
predicts a second-order phase transition, which indicates that ferromagnetic transition in dilute Fermi gases is 
possibly a counterexample to the BKV argument. The predicted critical gas parameter (k F a) c = 0.858 is in good 
agreement with the recent quantum Monte Carlo result (kpa) c = 0.86 for a nearly zero-range potential [S. Pilati, 
et al., Phys. Rev. Lett. 105, 030405 (2010)]. We also compare the spin susceptibility with the quantum Monte 
Carlo result and find good agreement. 

PACS numbers: 03.75.Ss, 05.30.Fk, 64.60.De, 67.85.-d 



I. INTRODUCTION 

Itinerant ferromagnetism is a fundamental problem in 
condensed-matter physics, which can be dated back to the ba- 
sic model proposed by Stoner . While the problem of itiner- 
ant ferromagnetism in electronic systems is quite complicated 
and the phase transition theory is still qualitative, a dilute spin- 
l Fermi gas with repulsive interactions may serve as a clean 
system to simulate the Stoner model. It is generally thought 
that the repulsive Fermi gas could undergo a ferromagnetic 
phase transition (FMPT) to a spin-polarized state with in- 
creased interaction strength H. Recently, the experimental- 
ists realized a two-component "repulsive" Fermi gas of 6 Li 
atoms in a harmonic trap by using a nonadiabatic field switch 
to the upper branch of a Feshbach resonance with a positive 
s- wave scattering length 00]. Therefore, it is possible to in- 
vestigate itinerant ferromagnetism in cold Fermi gases. The 
experimental progress in this direction has attracted intense 
theoretical interest JH-HH . 

The physical picture of the ferromagnetism in repulsive 
Fermi gases can be understood as a result of the competition 
between the repulsive interaction and the Pauli exclusion prin- 
ciple. The former tends to induce polarization and reduce the 
interaction energy, while the latter prefers balanced spin pop- 
ulations and hence a reduced kinetic energy. With increasing 
repulsion, the reduced interaction energy for a polarized state 
will overcome the gain in kinetic energy, and a FMPT should 
occur when the minimum of the energy landscape shifts to 



nonzero polarization or magnetization. 

Quantitatively, to study the FMPT in dilute Fermi gases at 
zero temperature, we should calculate the energy density & 
as a function of the spin polarization or magnetization x = 
(«f - «|)/(«-j- + «jj at given dimensionless gas parameter k?a 
which represents the interaction strength H ■ Here, fcp is the 
Fermi momentum related to the total density n = «|+«x by n = 
£p/(37r 2 ) and a > is the s-wave scattering length. Generally, 
the energy density can be expressed as 



8(x) = -nE F f(x), 



(1) 



where E-p = k 2 /(2M) is the Fermi energy with M being the 
fermion mass. The dimensionless function f(x), which de- 
pends on the gas parameter k F a, represents the energy land- 
scape with respect to the magnetization x. 

For the order of the FMPT, Belitz, Kirkpatrick, and Vojta 
(BKV) [3j] have argued that the phase transition in clean itin- 
erant ferromagnets is generically of first order at low tempera- 
tures, due to the correlation effects or the coupling of the order 
parameter to gapless modes that lead to a nonanalytic term in 
the free energy. The general form of the Ginzburg-Landau 
free energy for clean itinerant ferromagnets takes the form 



/glO) = ax 2 + ux 4 ln|jc| + fix 4 + 0(x 6 ), 



(2) 
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where we can keep /3 > 0. If the coefficient v is positive, the 
phase transition is always of first order. On the other hand, 
for negative v, one always has a second-order phase transi- 
tion. The BKV argument is based on the assumption v > 0, 
motivated by perturbation theory ||j7ll . This is true for many 
solid-state systems where the FMPT occurs at weak coupling. 
However, for dilute Fermi gases where the critical gas parame- 
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ter is expected to be of order 0(1), the assumption of a positive 
v is not reliable. 

In this paper, we study the nonperturbative effects on the 
FMPT by summing a set of particle-particle ladder diagrams 
to all orders in the gas parameter, motivated by the large- 
dimension expansion proposed by Steele ll38tl . We consider a 
universal repulsive Fermi gas where the effective range effect 
can be neglected, corresponding to a two-component "upper 
branch" Fermi gas with a positive s-wave scattering length. 
The prediction may also apply to the hard-sphere gas since 
the effective range corrections are sub-leading-order contribu- 
tions in the large-dimension expansion. Our main conclusions 
for the order and the critical gas parameter of the FMPT can 
be summarized as follows: 

(1) Order of phase transition. We predict a second-order 
phase transition, in contrast to the BKV argument. This sug- 
gests that the FMPT in dilute Fermi gas may correspond to the 
case of negative v. 

(2) Critical gas parameter. We predict a critical gas parameter 
(&f«)c = 0.858 where the spin susceptibility \ diverges. The 
critical gas parameter and the spin susceptibility we obtained 
are in good agreement with the quantum Monte Carlo results 

a. 

The paper is organized as follows. In Sec. II we briefly 
review the perturbative predictions for FMPT in dilute Fermi 
gases. In Sec. Ill we introduce the effective field theory ap- 
proach to the two-body scattering problem and show how we 
can recover the scattering amplitude by ladder resummation. 
We study the nonpertuabtive effects on FMPT in the theory 
of ladder resummation in Sec. IV and investigate the role of 
hole-hole ladders in Sec. V. We summarize in Sec. VI. 



II. PERTURBATIVE PREDICTIONS 

In the perturbation theory, the gas parameter k F a is treated 
as a small number. Up to the order 0((k F a) 2 ), the expression 
for f(x) is universal, that is, independent of the details of the 
short-range interaction. We have 



1 



fix) = -(n\ + ?il) + 



10k F a 
~9n~ 



3„3 



{ha)' 
2ln 2 



where r/i = (1 + x) 1/3 and r/i = (1 - x) lli . The zeroth-order 
term corresponds to the kinetic energy, and the first-order term 
coincides with the Hartree-Fock mean-field theory J2|. The 
coefficient ^(r]-\,m) in the second-order term was first evalu- 
ated by Kanno 13911 . Its explicit form is 

f = 22^(j7t + r]0 - 4i#n2L±i?i _ 4 ^ln-' 7T 



'7r 



(4) 



Setting x — 0, we recover the well-known equation of state for 
dilute Fermi gases, 



1 H kva 

9n 



4(11 -21n2) 
2~hfi 



(k F aY 



(5) 



which was first obtained by Huang, Yang, and Lee 14011 and 
recovered by Hammer and Furnstahl ll4lll in recent years using 
effective field theory. 

In the first-order perturbation, the FMPT is of second order 
and occurs at kpa = n/2 0]. However, taking into account 
the second-order corrections, one finds a first-order FMPT at 
kpa = 1.054 oi m. This can be understood by noticing 
the nonanalytical term oc jc 4 Ln| jc| with positive coefficient in 
the small-x expansion of the coefficient £(?7f , 77 j,). The small-x 
expansion of the energy density (3) takes the form 

f( x ) = /(0) + ax 2 + vx 4 \n\x\ + fix 4 + 0(x 6 ), (6) 

which is consistent with the BKV argument l37tl that the cor- 
relation effects or the coupling of the order parameter to gap- 
less modes generally leads to nonanalytical terms in the free 
energy. The coefficient u can be evaluated as 



40(k F a) 2 
243tt 2 



(7) 



Therefore, up to the order 0((k F a) 2 ), the Fermi gas problem 
corresponds to the case v > which is assumed in the BKV 
argument. 

In general, we expect that the critical parameter is of order 
0(1). Therefore, the perturbative predictions for the FMPT 
are probably unreliable. There naturally arises a serious prob- 
lem: Does the dilute Fermi gas problem really correspond to 
the case v > if the nonperturbative effects at k F a ~ 0(1) are 
taken into account? For the two-body problem in the vacuum, 
it is well known that an infinite set of bubble diagrams with the 
leading-order contact interaction must be resummed in order 
to reproduce the correct scattering amplitude if the two-body 
scattering length is large j42ll . Therefore it is natural to ex- 
tend the idea of ladder resummation to finite density so that 
the predicted equation of state works well even at k F a ~ 0(1). 
We can also compare the non-perturbative predictions with the 
results from recent quantum Monte Carlo (QMC) simulations 



m. LADDER RESUMMATION FOR TWO-BODY 
SCATTERING 

Before we establish a nonperturbative description for the 
FMPT in dilute Fermi gases, it is instructive to start with the 
low-energy effective field theory in vacuum ll4lll42Tl and to see 
how the two-body scattering amplitude is reproduced from the 
ladder resummation method. 

For nonrelativistic two-body scattering in the s-wave chan- 
nel associated with a short-range interaction, the scattering 
amplitude Jl(k) is related to the s-wave scattering phase shift 
6 by 



m.(k) 



An 



1 



M k cot 6 - ik 



(8) 
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where k = |k| is the scattering momentum in the center-of- 
mass frame. If there exist bound states for attractive interac- 
tions, the scattering amplitude should exhibit some imaginary 
poles, k = i -\J-ME\„ on the complex k plane with E\, < be- 
ing the binding energy. In general, the short-range interaction 
is characterized by a momentum scale A. Therefore, for low- 
energy scattering, that is, k <K A, the quantity k cot 5 can be 
expanded as a Taylor series in k 2 /A 2 . In quantum scattering 
theory, this is called the effective range expansion. 



kcotS 



i i °° / ; 2 V 

_I + Iy rnA 2 n 

1 1 ,2 

a 2 



i \«+i 



(9) 



where a is the scattering length and ro is the effective range. 
For a natural system II42I1 . we have \a\ ~ 1/A and |r„| ~ 1/A. 
An example commonly studied is a hard-sphere gas with ra- 
dius R, in which case a — R and ro = 2R/3. For cold atomic 
gases the interatomic interaction can be tuned by means of the 
Feshbach resonance, and we can have \a\ » |r„| ~ 1/A. 

According to the effective range expansion, one can con- 
struct the low-energy effective field theory ll4lll42ll describing 
scattering at momenta k <K A. Since we assume k <sz A, all in- 
teractions in the effective Lagrangian are contact interactions. 
The low-energy effective Lagrangian contains infinite contact 
interaction terms and is given by fiHl 



eft 
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id, + 



2M 

1 

(ifnlrf (<AVV) + H.c 



where Co and C2 are dimensionful coupling constants, V 

V Visa Galilei invariant derivative, and . . . denotes in- 
teractions with more derivatives (~ V 2 ", n > 2) which gen- 
erally have coupling constants C2,,. The coupling constants 
C2 n (n = 0, 1, 2, ...) should be determined by reproducing the 
scattering amplitude ZR&k). 

In practice, we can reproduce the scattering amplitude J[{k) 
order by order in a Taylor expansion in k/A. For small scat- 
tering length (\a\ ~ 1/A and \ak\ <K 1), we can expand the 
scattering amplitude as 



Ana 



1 — iak 



(ar 2 \ 
\-- a ) 



k 2 + . 



(11) 



However, for large scattering length (\a\ » 1/A), Kaplan et 
al. showed that one needs to expand Jl(k) in powers of k/A 
while retaining ak to all orders |42| : 



Mk) 



4n 



1 



M l/a + ik 



1 + 



ro/2 
1 / a + i k 



-k 2 + 



(12) 



This means if the scattering length is large, the loop diagrams 
with the leading-order interaction Co have to be resummed. 

According to the free fermion propagator £?o(/?o,p) = 
l/(po _ w p + ie) with the free dispersion a> p = p 2 /(2M), the 
(10) one-loop bubble diagram Bq(Pu, P) [FigQJ] is given by 



BoCPo,P) 



. f d 4 g 

J (2tt) 4 ^ + 



9o 



(P/2+q) 2 • Po 
2M T ,fc 2 



■<7o 



(p/2-q) 2 

2M 



+ ie 



r d\ 1 



(13) 



Here Po and P are the total energy and momentum of the pairs 
in the bubble diagram [see FigQla)]. Let p! and p 2 be the 
momenta of the scattering fermions, and we have P = pi + P2 
and k = (pi - p2)/2. Further, if the on-shell condition Po = 
(Pi + P?)/( 2M ) = P 2 /(4M) + k 2 /M is imposed, we find that 
Bo depends only on the relative momentum k, corresponding 
to the translational invariance. 

The integral over q in Eq. ( ITJb is linearly divergent and 
therefore needs to be regularized. A natural regularization 
scheme is to use a momentum cutoff equal to A 14311 . In this 
paper, we employ the dimensional regularization scheme. To 
this end, we change the space-time dimension from 4 to D and 
multiply the integral by a factor (/u/2) D . Here ji is an arbitrary 
mass scale introduced to allow the couplings C2,, multiplying 
operators containing V 2 " to have same dimensions for any D. 
In general, the integral Bo in D dimension can be evaluated as 



R(P P , _ r /3-D\ W2) 4 - D 

B„(p ,P) - -r^— j^zr^ 

/ p2 \(0-3)/2 

x M I -MPo + — -ie\ . (14) 

For small scattering length, it is convenient to use the minimal 
subtraction (MS) scheme which subtracts any l/(D - 4) pole 
before taking the D — > 4 limit. However, for large scatter- 
ing length, it is more convenient to use the power divergence 
subtraction (PDS) scheme. The PDS scheme involves sub- 
tracting from the dimensionally regularized loop integrals not 
only the 1 /(D - 4) poles corresponding to log divergences, as 
in MS, but also poles in lower dimensions which correspond 
to power law divergences at D = 4. The integral Bo has a pole 
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in D = 3 dimensions. It can be removed by adding a countert- 
erm 6B = Mfi/[An(3 - D)] to B El- Finally, the subtracted 
integral in D — 4 dimensions is 



B (P ,V) = - — 



M 
An 



p2 

-MP + — 



(15) 



Note that the MS scheme corresponds to the p — case. 

The dependence of C2„(yu) on p is determined by the re- 
quirement that the scattering amplitude is independent of the 
arbitrary mass scale p. To this end, we impose the on-shell 
condition, P = P 2 /(4M) + k 2 /M. Then the one-loop bubble 
diagram becomes Bo(k) = -M(p + ik)/(An). Summing the 
bubble diagrams with Co vertices, we obtain ll42tl 

C 2 (v)k 2 



Co(p) 



1 - C (fi)B (k) [i _ C Q (p)B Q (k)] 2 
Comparing this result with the expansion (fTZt . we obtain 
An 1 



(16) 



C (m) = 



C 2(M) = 77 



M -fi + l/a' 
An I 1 
M 1/a 



2 ' 



(17) 



It was shown that these results fulfill the renormalization 
group equations ll42tl . We note that the mass scale p is sim- 
ilar to the cutoff A. In the cutoff scheme, we have Co(A) = 
(47r/M)(-2A/7r + 1/a) 0. 

In the following, we mainly consider a short-range poten- 
tial with a positive scattering length a and negligible effective 
range ro <s a. In this case, we are able to obtain a univer- 
sal result for fix) which is independent of the details of the 
interaction. In this case, the pair propagator Sq(Pq, P) in the 
vacuum is given by 

CoOu) 



1 - Co(M)B (Po, P) 
An 1 



M 



l/a- tJ-MP + f - 



(18) 



For positive scattering length, the pair propagator has a pole 
given by Pq = -l/(Ma 2 ) + P 2 /(4M). This pole corresponds 
to a bound state with binding energy E h = -1 1 (Ma 2 ) and 
effective mass 2M. Therefore, if the effective range is neg- 
ligible, the underlying potential must be attractive and the 
ground state is a bound molecule of size a. However, for 
two-body scattering state with positive center-of-mass energy 
E = k 2 /M > 0, the effective force is repulsive. This is the so- 
called "upper branch," which is well defined in the two-body 
picture. For the many-body problem, a metastable "repulsive" 
Fermi gas can be realized if all fermions are forced on the up- 
per branch of a Feshbach resonance with a positive s-wave 
scattering length flHl. 



IV. LADDER RESUMMATION AT FINITE DENSITY 

We now turn to the many-body problem of the repulsive 
Fermi gases. The main purpose of this paper is to present a 



(a) 



(b) 




Pi 



B 



FIG. 1: (a) The elementary particle-particle bubble B(p u p 2 ) with 
external momenta p\ and p 2 for the two spin components. The solid 
line with arrow corresponds to the particle term of the propagator 
H9\ . The dashed line represents the interaction vertex Co- (b) A 
typical particle-particle ladder diagram contributing to the interaction 
energy. 



nonperturbative calculation for function f(x) from which we 
can give a better prediction for the FMPT In general, we ex- 
pect that the new nonperturbative result for f(x) satisfies the 
following two criteria: (i) The function f(x) recovers the per- 
turbative result, Eq. (0, at weak coupling kpa — > 0; and (ii) 
since we consider a short-range potential with ro <K a, the 
physical result should be universal, i.e., f(x) depends only on 
the gas parameter k^a and does not depend on other parame- 
ters such as the renormalization scale //. The criterion (ii) is 
hard to fulfill since the loop corrections in quantum field the- 
ory generally bring the renormalization scale dependence and 
weaken the prediction power due to the arbitrariness in the 
choice of the renormalization scale. However, in the follow- 
ing we show that the result from the particle-particle ladder 
resummation, corresponding to the leading order of the large- 
dimension expansion, is independent of the arbitrary mass 
scale pi. 



A. Nonperturbative energy density 

We first construct the nonperturbative version of the energy 
density &(x) using the vertex Co(p) determined in Sec. Ill 
and the free propagators for the two spin components at finite 
density [45] 

Q(|p|-*g) 0(*g-|p|) 

0Apo,p)= V + -. (19) 

po - ojp + le po - o) p - ie 

Here cr =f , J,, fep = kprj-^ i are the Fermi momenta of the two 
spin components and 0(z) is the Heaviside step function. For 
each spin component, the propagator (fT9] l describes two types 
of excitations, particles with momentum |p| > k^ and holes 
with |p| < k%. 

The dilute imperfect Fermi gases are best described by re- 
summing the multiple interactions in terms of the scatter- 
ing amplitude. The Galitskii integral equations [45] for the 
effective two-particle interaction or scattering amplitude in 
the medium are given by the ladder resummation ll46ll . On 
the other hand, for large gas parameter kpa, one may look 
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for other expansion parameters instead of k-pa itself. Steele 
ll38tl and Schafer et al. j47[| have suggested a new expan- 
sion method using 1/2) as the expansion parameter, where 
D = 2 D I 2 with D being the space-time dimension. Most im- 
portant, they have shown that the contribution of the particle- 
particle (pp) ladder resummation, fiff , i s the leading-order 
contribution of the 1/2) expansion yl|,|47|], that is, 

6 = 6 Bn +^ + 0(l/2>). (20) 

All other contributions like hole-hole (hh) ladder sum and ef- 
fective range corrections are suppressed by a factor 1 /2). Ac- 



cording to the above arguments, we expect that the most im- 
portant nonperturbative contributions come from the leading 
order of the 1/2) expansion. The interaction energy density 
£^ to this order is given by all the particle-particle scattering 
terms (i.e., the npp-lhh bubbles for all n — 0, 1, 2, • ■ ■ ). 

To evaluate the interaction energy density, we first calculate 
the elementary in-medium particle-particle bubble B(Po,P) 
shown in Fig. [TJa). The fermion lines in the bubble dia- 
gram correspond to the particle terms of the free propagator 
(fT~9b . According to the finite-density Feynmann rules B41I1 . it 
is given by 



B(P ,P) 



. r d*g W 
J (2tt) 4 ^ + 



0(|P/2 + q| 



0(|P/2-q|-4) 



(P/2+q) 2 , ,• Po 

d\ @(|P/2 + q| - 4)0(|P/2 - q| - k l F ) 
{2nf 



•s 



(P/2-q) 2 
2M 



- 2- + ie 

AM M 



(21) 



For vanishing densities, k F = 0, the in-medium particle-particle bubble recovers the vacuum result Bo- If the on-shell condition 
is imposed, the in-medium particle-particle bubble B depends on not only the relative momentum p but also the total momentum 
P. This is due to the loss of translational invariance in the presence of Fermi sea. 
We can separate B into a vacuum part and a medium part using the identity 

0(|P/2 + q| - 4)©(|P/2 - q| - 4) = 1 - ®(k\ ~ |P/2 + q|) - 0(4 - |P/2 - q|) + 0(4 - |P/2 + q|)0(4 - |P/2 - q|). (22) 

The vacuum part (corresponding to 1) is identical to Bo defined in the last section and is linearly divergent. The medium part 
is convergent. For the vacuum part, it is natural to use the dimensional regularization with PDS scheme introduced in the last 
section. 

Then the zipp-lhh bubble [see Fig. [TJb) for a typical example] at given n reads 



) 4 En 



0(4 - |P/2 + k|) 0(4 - |P/2 - k|) 



, (P/2+k) 2 • P 
+ k ° 2M T 



_ ^0 2M le 



[B(Po,P)T 



f _r-i f d4p f dA> 



r 



where e"' p ° with r/ — > + is a convergence factor 114 ill . The 
integration over Po picks up the pole or imposes the on-shell 
condition P = P 2 /(4M) + k 2 /M, that is, 



dP j nPo [B(P ,V)r 
2ni 6 Po - — 



M Lt 



[B(p,k)r, 



where the on-shell version of B is given by 
B(P, k) = 

d\ 0(|P/2 + q| - 4)©(|P/2 - q| - k l F ) 



M 



J (2nf 



k 2 - q 2 + ie 



The total interaction energy density 6[? is given by 



n=0 



(24) 



(25) 



(26) 



Completing the summation of this geometric series, we obtain 
the interaction energy density at the leading order of the 1/2) 
expansion, 



o(0) = c f d Pi f d & 
int °J (2tt) 3 J {2nf 



d 3 p 2 0(4 -|pi 1)0(4 -|p 2 |) 



l-Cofi(P.k) 



(27) 



where pi,2 = P/2 ± k as defined in Sec. III. The imaginary 
part of B can be evaluated as 



ImB(P,k) = -^Me(|pi| -4)0(|p 2 | -4). 



(28) 



This quantity is nonzero only when the momenta pi and p 2 
are both above the Fermi surfaces. However, the final inte- 
gration over p! and p 2 in the interacting energy density £i nt is 
associated with a phase-space factor 0(4 ~ IPi 1)0(4 _ Ip2l)- 
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Therefore, the interaction energy density is real and physical, 
as we expected. 




FIG. 2: The pole energy E (divided by the binding energy |£t,l = 
1 1 {Ma 1 ) in the vacuum) at zero pair momentum P = as a function 
of the gas parameter kpa. The pole energy turns out to be positive for 
kpa > n/2. 




FIG. 3: (Color online) The pole energy E (divided by 2Ep) as a func- 
tion of the pair momentum P = |P| (divided by kp) for various values 
of the gas parameter kpa. The dashed line corresponds to the disper- 
sion E(P) = P 2 /(4M). 



B. In-medium two-body problem 

Since we adopt a zero-range potential, a bound state with 
binding energy E\> = -1 /(Ma 2 ) always exists for positive 
scattering length a > 0. A key problem here is that how 
we can describe a metastable repulsive Fermi gas where all 
fermions are forced on the scattering states. Actually, the so- 
called upper branch has clear meaning only in the two-body 
picture, and so far it is not clear to what extent this two-body 
picture of a "repulsive" Fermi gas will persist. A recent study 



of three attractive fermions shows that there are many non- 
trivial avoided crossings between the two branches close to 
the resonance (a — > oo), making it difficult to unambiguously 
identify a repulsive Fermi system lfl7ll . 

To realize a metastable repulsive Fermi gas we have to ex- 
clude the molecule bound states of two atoms with unlike 
spins and enforce all atoms to the scattering states 04811 . One 
possible prescription is to subtract the contribution from the 
bound-state poles within the Nozieres-Schmitt-Rink (NSR) 
theory ll48tl . However, as designed, NSR theory works well 
only at temperature higher than the critical temperature of 
superfluidity. The particle-particle resummation theory we 
present complements the NSR theory and can be regarded as 
the zero-temperature analog of the NSR theory. 

The key point in this problem is to consider the medium ef- 
fects on the bound-state properties. To this end, we first con- 
struct the pair propagator S(Pq, P) in the presence of Fermi 
seas. With the in-medium elementary particle-particle bubble 
B(Pq, P), the in-medium pair propagator S(Po, P) is given by 
the ladder resummation, 



S(P ,F) 



1 - C Qi)B(P ,P) 



(29) 



With this pair propagator, the interaction energy density £ 
can be expressed as 



(0) 



? (0) 



r d s ~p r 



IPll)®(4-|P2l) 



J 2^ 



S(P ,P) 



"0 — 777 — 77 



(30) 



In general, the in-medium pair propagator S(Pq, P) has a 
real pole Pq = E(P) corresponding to the in-medium bound 
state. However, we now show that such pole does not con- 
tribute to the energy density for the regime of the gas param- 
eter k F a we are interested in. As shown in the last subsection, 
in the calculation of the interaction energy density £j nt , the on- 
shell condition Pq = P 2 /(4M) + k 2 /M is imposed and the in- 
tegrations over the momenta P and k are performed according 
to the finite-density Feynmann rules. Therefore, if the energy 
dispersion of the pole E(P) satisfies the condition 



p2 

E(P) < 

V ; 4M 



(31) 



for arbitrary P, its contribution to the energy density is natu- 
rally excluded. 

Since the main purpose of this paper is to study the 
FMPT which corresponds to an instability toward a small 
polarization x, we can set x = here. For convenience, 
we define two dimensionless quantities s = \P\/(2kp) and 
z = y/MP - P 2 /4 + ie/kp = ^P /(2E F ) - s 2 + ie. The in- 
medium pair propagator can be evaluated as 



S(P ,P) 



An 



1 



M l/a-(k F /jT)W(s,z) 



(32) 



where W(s, z) is given by 
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W(s,z) = 



1+s-z 1 
1 + s + z In h 



1 + s + z 



s 2 -z\ (l + s) 2 -z 
In 



2s 



(l-z) 2 -s 2 l-/-z 2 , 
2 + z In — — — + = In 



1 - s 2 - z 2 
(l+s) 2 -z 2 



s 2 -(l +z) 2 



2s 



(l-s) 2 -z 2 



©(1 - s) 

0(5-1). 



(33) 



The same result was also obtained in a recent paper ll49ll . 

For zero pair momentum P = 0, the condition £(P) < 
P 2 /(4M) implies £(0) < 0. We thus focus on the regime of 
the gas parameter where the pole £(0) is negative. In this case, 
£(0) is determined by a simple equation, 

7T, — = 1 + a/ -TE— arctan V ~TE~- (34) 
2k F a V 2£ F V 2£ F 

This equation has negative solution only for < kpa < tt/2, 
where the solution represents the binding energy of a in- 
medium bound state. The numerical result for £(0) is shown 
in Fig. [2] In the low-density limit kp — » 0, £(0) recovers 
the vacuum result £ vac (0) = £\, = -I /(Ma 2 ). However, at 
finite kp, the medium shields the bound state and reduces the 
binding energy, that is, |£(0)| < |£'vac(0)|. For kpa > n/2 and 
kpa < 0, Eq. (f34-b has a positive solution which corresponds 
to the positive energy pole of the in-medium pair propagator. 
Such a pole is associated with Cooper pairs, and its appear- 
ance represents the BCS instability. This positive energy pole 
does not lead to singularities in the energy density integration, 
as can be seen in next subsection, and does not need special 
treatment. 

Note that £(0) < is not a sufficient condition for £(P) < 



P 2 /(4M). We thus have to check the energy dispersion £(P) 
carefully. The numerical results for some values of the gas pa- 
rameter kpa are shown in Fig. [3] For kpa < 1 .34, the condition 
£(P) < P 2 /(4M) is fulfilled for all values of P. However, for 
1.34 < kpa < n/2, there exists a regime P\ < |P| < £2 where 
£(P) > P 2 /(4M). 

In conclusion, the condition £(P) < P 2 /(4M) is fulfilled for 
kpa < 1.34. Therefore, in the parameter regime kpa < 1 inves- 
tigated in the following, the contribution from the bound state 
can be naturally excluded in the ladder resummation scheme. 



C. Evaluating the energy density 

Now we evaluate the explicit form of the energy density 
£(x) and the dimensionless function f(x). First, the elemen- 
tary particle-particle bubble B(P, k) can be decomposed into 
four parts 

B(P, k) = fl (P, k) + £ T (P, k) + B X (P, k) + %(P, k), (35) 

where Z?o is the vacuum part discussed in Sec. Ill and the other 
parts are given by 



C d'q &(k' - \P/2 + q|) r d\ ®(kl - |P/2 - q|) 

B T (P,k) = -M - * \ F 2 ' 2 . 4l; , B i (P,k) = -M - ^ 3 \ F 9 [ 
J (2ny k 2 - q 2 + le J (2ny k 2 - q 2 + 

B n (P,k) = M — — J— ■ 

J (2jiy k 2 -q 2 + ;e 



(36) 



For convenience, we define another dimensionless quantity t = |k|/fcp together with s as defined in the last subsection. Since 
the imaginary part of B does not contribute to the interaction energy, we need only to evaluate the real part of B. We have 

ReB (M) = -^. Refl T (s,0 = ^f*T(J.O, ReB A (*,/) = ^f*±(M), ReB n (s,t) = ^-R n (s,t). (37) 



where Ro-(s, f) (cr =f , J,) reads 

R a (s,t) = 
and the function Rfi(s, t) is 

Rn{s,t) = \ 



■ (s + ty 



4s 



In 



rjo- + s + t 



TJa- - S-t 



ril-(s-tf 
4s 



In 



TJtr + S-t 



s + t 



+ llo-, 



(38) 



-Q(x)R l (s, t) - ®(-x)Rj(s, t) , < s < i|77 T - 77J 

K 1 (s,t) + K L (s,t) , i|77 T -77J < s < i|j/ T +774,1 
, elsewhere. 



(39) 
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Here A"o-(s, t) is defined as 



As 



-In 



for - ^) 2 - f 2 



r 2 _ ,2 _ ,2 



+ 2 ln 



77o- - 5 + t 



S-TJ 



- S-t 



(40) 



where r = (77^ + ?7 2 )/2. 

Finally, the elementary particle-particle bubble reads 

Mu Mkp 

B(s,t) = + -^-R w {s,f), (41) 

where the function R pp (s, t) is defined as 

R pp (s, t) = R t (j, f) + f) + « tt (j, 0, (42) 

Substituting the result of B(s, t) into the expression of we 
observe that the energy density is independent of the renor- 
malization mass scale fi. Converting the integration variables 
Pi and P2 to P and k, we find that the function f(x) can be 
expressed as 

1 so f °° r°° 

fix) = ^ + — J * 2 ds J Q tdtI(s,t)F(s,t), (43) 
where F(s, t) is given by 
F(s,t) = 



kpa 



1 - j-k F aR pp (s, t) 



(44) 



The function I(s, t) appears due to integration over the angle 
between P and k. Its explicit form is 



I(s,t) = 



tf-(s + t) 1 



As 



-®(s + t-7]<t) + (t7 T -> 774.) + t 
f 2 )0(/7 T -| S -f|)0(^-k-f|). (45) 



As we mentioned in the beginning of this section, it is im- 
portant to check whether the present result for f(x) is con- 
sistent with the perturbative expression (|3) for weak coupling 
kpa <k 1. To this end, we expand the function F(s, k) as 

F(s, t) = k F a + -(k F a) 2 R pp (s, t) + 0((k F a) 3 ). (46) 

71 

Using the expressions for I(s, f) and R pp (s, t), we can show 
that 

so r°° , r°° 10 , , 

— s z ds tdtl(s, t) = — 77^77; (47) 
n Jo Jo 9?r 1 



and 



so p 2 r 



tdtl(s, t)R pp (s, t) = 



2ln 2 ' 



(48) 



Therefore, our nonperturbative expression (27) exactly recov- 
ers the perturbative result (3) at weak coupling. This con- 
vinces us that the present theoretical approach is suitable to 
study the universal upper-branch Fermi gas with a positive 
scattering length. In addition, we can compare the results from 
our theory and the second-order perturbation on the same foot- 
ing and study the nonperturbative effects on the FMPT 



D. Results and discussion 

(A)Energy density and compressibility. We first study the 
equation of state for the unpolarized case x = 0. The gas 
parameter dependence of the energy density & in the regime 
< k F a < 1 is shown in Fig. [4] We find that the result from 
the ladder resummation is consistent with the perturbative re- 
sult (3) for small gas parameters k F a < 0.4. However, signifi- 
cant deviations are found for k F a > 0.4, consistent with recent 
quantum Monte Carlo simulations JH @]. For the quantum 
Monte Carlo simulations of the attractive interactions with a 
negligible effective range (corresponding to UB and UB2 in 
Fig. |4|i, the exclusion of molecular bound states is imple- 
mented by choosing a two-body Jastrow factor iH |^] to be 
the scattering solution of the attractive potential correspond- 
ing to positive energy, which, by construction, is orthogonal 
to the bound molecules. Therefore, the accuracy of the quan- 
tum Monte Carlo data depends on the choice of the Jastrow 
factor. Actually, exact orthogonality of the many-body vari- 
ational wave function to the superfluid ground state (molecu- 
lar condensation) can not be achieved in the quantum Monte 
Carlo simulations We note that our theoretical curve 

agrees better with the UB2 data than with the UB data. The 
reason could be that the UB2 data from [6J are obtained with 
a Jastrow factor which imposes a better orthogonality to the 
superfluid ground state. 

An important issue is whether the system is mechanically 
stable. The mechanical stability of the system requires a pos- 
itive compressibility k, which is defined as 



1 2 d 2 & 
k on- 



(49) 



For the present ladder resummation theory, the explicit form 
of k can be evaluated as 

— = 1 + I s 2 ds tdtl(s, t)G(s, t), (50) 

« n Jo Jo 

where kq = 3/(2nE F ) is the compressibility for noninteracting 
Fermi gases, and the function G{s, t) is given by 



G(s, t) = F(s, t) + 



5R pp (s,t)^ 2/ _ flpp(M)^ 



9n 



-F z (s,t) + 



9jt 2 



F 3 (s,t). (51) 



The compressibility k as a function of the gas parameter k F a 
is shown in Fig. [5] Comparing to the result from the second- 
order perturbation theory, 

ko ,2, 8(ll-21n2), f . 2 

— = 1 + -k ¥ a + — (kpaf, (52) 

k n 157r z 

good agreement is found for small gas parameters, as we ex- 
pected. In the regime < kpa < 1 we are interested in, we 
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FIG. 4: (Color online) The energy density £ (divided by its value 
Go = InZsp for noninteracting Fermi gases) as a function of the gas 
parameter kpa (0 < kpa < 1) for the unpolarized case x = 0. The 
solid line is the result calculated from our particle-particle ladder re- 
summation theory (RTh). The dashed line is result of the second- 
order perturbation theory (PTh). The dash-dotted horizontal line 
corresponds to the energy of the fully polarized state (x =1), i.e., 
f(l) = 2 2,i . The blue squares are the quantum Monte Carlo (QMC) 
data for the hard sphere (HS) potential Q, the red circles are for the 
upper branch (UB) of a square well potential @], and the green di- 
amonds are for the upper branch (UB2) of an attractive short range 
potential @]. For UB and UB2 cases, the effective range r is much 
smaller than the s- wave scattering length a |@,@]. 



"2 




FIG. 5: The compressibility k [divided by its value ko = 3/(2nEp] 
for non-interacting Fermi gases) as a function of the gas parameter 
kpa (0 < kpa < 1) for the unpolarized case x = 0. The solid line 
is the result calculated from our particle-particle ladder resummation 
theory. The dashed line is result of the second-order perturbation 
theory. 



find that the compressibility k is positive, indicating that the 
system is mechanically stable. 

(B)Spin susceptibility. Next we study the response of the 



energy density to an infinitesimal polarization x. This re- 
sponse is referred to as the spin (or magnetic) susceptibility. 
The spin susceptibility x can be defined as 



I - 1 d 2 £ 

X n 2 dx 2 



3E F d 2 f(x) 



x=0 



5n dx 2 



(53) 



A=0 



In the present ladder resummation theory, an explicit form of 
X is hard to obtain. In practice, we expand the function f(x) 
near x = as f(x) = f(0) + ax 2 + ■ ■ ■ . The coefficient a is 
related to the spin susceptibility by 



Xo 
X 



(54) 



where xo = 3« / (2Ep) is the spin susceptibility of noninter- 
acting Fermi gases. Therefore, a diverging spin susceptibility 
generally indicates a FMPT, as long as the transition is of sec- 
ond order. 

In the second-order perturbation theory, an analytical result 
for^- can be achieved, 



x 



l 



> 16(2 + ln2) 

KpU \K-F U ) I 

r 15n~ 



(55) 



which indicates a diverging spin susceptibility at kpa = 1.058. 
However, this differs from the critical gas parameter {kpa) c = 
1.054, because the phase transition is of first order in the 
second-order perturbation theory due to the appearance of the 
nonanalytical term vx A In |jc| with v > 0. 

Our result for the spin susceptibility x as a function of the 
gas parameter kpa is shown in Fig. [6] and compared with the 
perturbative result. We find that the spin susceptibility pre- 
dicted by the ladder resummation deviates significantly from 
the second-order perturbative result for kpa > 0.4. Further, 
the spin susceptibility diverges at kpa = 0.858, in contrast 
to the value 1.058 from the second-order perturbation theory. 
The data from the quantum Monte Carlo simulations JH are 
also shown in Fig. [6] as a comparison. Our theoretical re- 
sult is in good agreement with the data for the upper branch 
of the square well potential where the effective ranee ro is 
tuned to be much smaller than the scattering length H. The 
gas parameter kpa = 0.86 where x diverges is very close to 
our prediction kpa = 0.858. For the purely repulsive po- 
tential, that is, the hard-sphere potential, the effective range 
effect cannot be neglected a priori. However, we find that 
our result still has nice agreement with the data for the hard- 
sphere case. The gas parameter kpa = 0.82 where^ diverges is 
also close to our prediction kpa = 0.858. Actually, the differ- 
ence between the upper branch and the hard sphere cases [i.e., 
0.86 - 0.82 = 0.04] is very small compared with the critical 
gas parameters. This indicates that the contribution from the 
effective range effect is relatively small even for kpa ~ 0(1), 
if the quantum Monte Carlo results are reliable. This can be 
understood from the large-dimension expansion ll3~8l l47ll in- 
troduced in the beginning of this section: The particle-particle 
ladder sum is the leading-order contribution in the 1 /£) expan- 
sion, and all other contributions including the effective range 
corrections are suppressed by a factor 1 I'D. 
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der 0((kpa) 3 ) even for the unpolarized case x = in the pertur- 
bation theory iRHl . However, some definite conclusions can be 
drawn from our numerical results: (1) Higher-order terms in 
the gas parameter can also generate nonanalytical terms of the 
form x 4 ln|x| and may generate other important non-analytical 
terms which are not known due to the mathematical limita- 
tion. (2) The coefficients of the nonanalytical terms generated 
by the higher-order contributions are certainly not always pos- 
itive, and they are generally proportional to (kpa)" for the nth- 
order contributions. Since the phase transition occurs at a gas 
parameter kpa ~ 0(1), the nonperturbative effects from the 
sum of the higher order contributions are very important. As 
we have shown numerically, their effects are not only reduc- 
ing the critical value of the gas parameter but also changing 
the order of the phase transition. 



k F a 



FIG. 6: (Color online) The dimensionless inverse spin susceptibility 
Xolx as a function of the gas parameter k f a. The blue squares and red 
circles are the QMC data H 

for the HS and UB cases, respectively. 
The solid line is the result calculated from the particle-particle ladder 
resummation. The dashed line is the second-order perturbative result. 



(C)Ferromagnetic transition. While a diverging spin sus- 
ceptibility indicates a ferromagnetic phase transition, the or- 
der of the ferromagnetic phase transition and the critical gas 
parameter (kpa) c should be obtained by studying carefully the 
shape of the energy landscape, that is, the full x dependence of 
the function f(x). To very high numerical accuracy, we have 
not found any maximum at x + in the energy landscape. In- 
stead, we find a second-order phase transition at kpa = 0.858, 
where the function fix) starts to develop a minimum at x + 0, 
consistent with the gas parameter where the spin susceptibil- 
ity diverges. This is in contrast to the second-order pertur- 
bation theory which predicts a first-order phase transition at 
k^a = 1.054 iffill . where the spin polarization a jumps from 
zero to x c = 0.573. A second-order FMPT for a zero-range 
potential model was also obtained by Heiselberg lfl8ll recently 
using a completely different many-body method. 

It seems that our result of a second-or der p hase transition 
is in contradiction to the BKV argument IB7I1 . However, the 
BKV argument is based on the assumption that v > 0. Ac- 
tually, we have fitted the energy density of the form f(x) = 
f(Q) + ax 2 + vxHn\x\ + fix 4 for small x. For small gas param- 
eter kpa < 0.3, the coefficient v agrees well with the pertur- 
bative result v = 4Q(kfa) 2 /(24-3n 2 ). However, for larger k^a 
(especially around the critical gas parameter), it turns out to 
be negative due to the nonperturbative effects. This indicates 
that the FMPT in the systems of dilute repulsive Fermi gases 
corresponds to the case u < and is a counterexample to the 
BKV argument where the assumption v > is adopted. 

Since an analytical expression for the function f(x) as well 
as the coefficient v cannot be achieved in the present ladder re- 
summation theory, we cannot understand analytically how the 
nonparturbative effects modify the order of the phase transi- 
tion. In fact, analytical results cannot be obtained from the or- 



FIG. 7: The elementary bubbles organized in the number of the ML 
The solid line with a cut represents the MI part of the propagator, and 
the pure solid line corresponds to the vacuum part. 



V. INCLUSION OF HOLE-HOLE LADDERS 

In this section we check whether our conclusion that the 
FMPT is of second order is changed by other contributions. 
We consider the contributions from the hole-hole ladder dia- 
grams by summing the combined particle-particle and hole- 
hole ladders to all orders in kpa while keeping the criteria (i) 
and (ii) satisfied. Such a resummation scheme for the unpo- 
larized case x — has been performed b y K aiser f50ll . 

Following the treatment by Kaiser 15011 . we rewrite the 
propagator ( fT~9b in an alternative form 

GApo, P) = QoiPo, P) + 2mS(p - w p )®(&£ - |p|), (56) 

where the first term corresponds to the vacuum propagator 
@o(po,p) = (po - w p + ie)' 1 and the second term is a so- 
called medium insertion (MI) ll50ll . The elementary bubbles 
in this treatment are shown in Fig. [7] The first diagram Bq 
is identical to the vacuum part studied in Sec. Ill and it can 
be renormalized using the PDS scheme. For our purpose of 
resummation, we are interested in the following two quanti- 
ties, Bo + B\ + Z?2 and Z?o + B\, which are mutually complex 
conjugate. We have 

Mu Mk F 

B + B l+ B 2 = --Z + —f[R(s,t)-i7rI(s,t)], 

Bo + Bi = -^tL + ^L [ R ( S , t) + frr/( Jt f)] t (57 ) 
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where R(s, t) = R-[(s, t) + Ry(s, t) and I(s, t) is the function 
defined in (30). 

To sum all ladder diagrams built from the elementary bub- 
bles, we first notice that the nonvanishing contributions to the 
interaction energy come from diagrams with at least two ad- 
jacent Mis f50ll . Then a typical nth-order contribution would 
look like the ring diagram of Fig. [U(b) with n vertices and at 
least two adjacent Mis. Naively, all these nth-order diagrams 
are summed to give g"[(B Q + Bi +B2T -{Bq + B\)"], where the 
subtraction gets rid of those diagrams which have no adjacent 
MI pairs. However, this expression is complex and therefore 
cannot be the correct one. The crucial observations are that as 
follows: (1) Each nth-order ring diagram has an n-rotational 
symmetry. Therefore, we should introduce an additional fac- 
tor 1/n. (2) An nth-order ring diagram comes from closing 
two open MI lines of an nth-order ladder diagram, which in- 
troduces an integration over the allowed phase space |pi| < k F 
and IP2I < k F but does not contribute a factor B2 to the energy 
as the naive expression does. These amendments lead to the 
correct nth-order contribution to the interaction energy ll50l : 
g"[(B Q + B x + B 2 f - (B + B{) n ]/(2ilri). The summation over 
n leads to two complex-conjugated logarithms and the final 
result is real. 

The final result for the energy density does not depend on 
the renormalization scale /j., and the function f(x) in this re- 
summation scheme also takes the form (l43l , while the func- 
tion F(s, t) becomes 

In [l - -k F aR(s, t) + ik F aI(s, t)] - c.c. 

F(s, f) = -t n - — — . (58) 

2il(s, t) 

For small gas parameter k F a <k 1, F(s, t) can be expanded as 

1 , , 
F(s, t) = k F a + -(k F a) 2 R(s, t) + 0((k F aY). (59) 
n 



We can also check that 



which reflects the fact that the hole-hole ladders start to con- 
tribute at the order 0((k F a) 3 ) JM EH El • Therefore, the cri- 
teria (i) and (ii) are also fulfilled in the present resummation 
theory. Numerically, we also find a second-order phase tran- 
sition, which occurs at a smaller gas parameter k F a = 0.786. 
We note that the inclusion of hole-hole ladders may not im- 
prove the quantitative result, since it only includes part of the 
beyond-leading-order contribution in the large-£) expansion. 



VI. SUMMARY 

In summary, we have studied the nonperturbative effects on 
the ferromagnetic phase transition in repulsive Fermi gases by 
summing the ladder diagrams to all orders in the gas parame- 
ter k F a. The nonperturbative effects not only reduce the criti- 
cal gas parameter but also change the order of the phase tran- 
sition. The resummation of particle-particle ladders, which 
corresponds to the leading order of the large-dimension ex- 
pansion, predicts a second-order phase transition occurring at 
k F a = 0.858, in good agreement with the quantum Monte 
Carlo result [5J. The spin susceptibility calculated from our 
resummation theory is also in good agreement with the quan- 
tum Monte Carlo results. Therefore, the resummation of the 
ladder diagrams provides a more quantitative way to study the 
ferromagnetic transition in repulsive Fermi gases. In this pa- 
per, we have considered only a zero-range potential model. It 
will be interesting to study the nonuniversal shape-dependent 
contributions using the finite-density effective range expan- 
sion |E3l. 



7T Jo Jo 



tdtl(s, t)R(s, t) = f^Ii^l , (60) 
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It is generally believed that a dilute spin- 1/2 Fermi gas with repulsive interactions can undergo a ferromagnetic 
phase transition to a spin-polarized state at a critical gas parameter (kpa) c . Previous theoretical predictions of 
the ferromagnetic phase transition are based on the perturbation theory which treats the gas parameter as a small 
number. On the other hand, Belitz, Kirkpatrick, and Vojta (BKV) have argued that the phase transition in clean 
itinerant ferromagnets is generically of first order at low temperatures, due to the correlation effects that lead to 
a nonanalytic term in the free energy. The second-order perturbation theory predicts a first order phase transition 
at (k F a) c = 1.054, consistent with the BKV argument. However, since the critical gas parameter is expected to 
be of order 0(1), perturbative predictions may be unreliable. In this paper we study the non-perturbative effects 
on the ferromagnetic phase transition by summing the particle-particle ladder diagrams to all orders in the gas 
parameter. We consider a universal repulsive Fermi gas where the effective range effects can be neglected, 
which can be realized in a two-component Fermi gas of 6 Li atoms by using a nonadiabatic field switch to the 
upper branch of a Feshbach resonance with a positive s-wave scattering length. Our theory predicts a second 
order phase transition, which indicates that ferromagnetic transition in dilute Fermi gases is possibly a counter 
example of the BKV argument. The predicted critical gas parameter (kpa) c = 0.858 is in good agreement with 
recent Quantum Monte Carlo result (kpa) c = 0.86 for a nearly zero-range potential [S. Pilati, et ah, Phys. Rev. 
Lett. 105, 030405 (2010)]. We also compare the spin susceptibility with the Quantum Monte Carlo result and 
find good agreement. 

PACS numbers: 05.30.Fk, 03.75.Ss, 64.60.De, 71.10.Ca 



I. INTRODUCTION 

Itinerant ferromagnetism is a fundamental problem in con- 
densed matter physics, which can be dated back to the basic 
model proposed by Stoner [ 1 ] . While the problem of itinerant 
ferromagnetism in electronic systems is quite complicated and 
the phase transition theory is still qualitative, a dilute spin- 1/2 
Fermi gas with repulsive interactions may serve as a clean sys- 
tem to simulate the Stoner model. It is generally thought that 
the repulsive Fermi gas could undergo a ferromagnetic phase 
transition (FMPT) to a spin-polarized state with increased in- 
teraction strength H. Recently, the experimentalists realized 
a two-component "repulsive" Fermi gas of 6 Li atoms in a har- 
monic trap by using a nonadiabatic field switch to the upper 
branch of a Feshbach resonance with a positive s-wave scat- 
tering length Jit]. Therefore, it is possible to investigate itin- 
erant ferromagnetism in cold Fermi gases. The experimental 
progress in this direction has attracted intense theoretical in- 
terest HUH. 

The physical picture of the ferromagnetism in repulsive 
Fermi gases can be understood as a result of the competition 
between the repulsive interaction and the Pauli exclusion prin- 
ciple. The former tends to induce polarization and reduce the 
interaction energy, while the latter prefers balanced spin pop- 
ulations and hence a reduced kinetic energy. With increasing 
repulsion, the reduced interaction energy for a polarized state 
will overcome the gain in kinetic energy, and a FMPT should 
occur when the minimum of the energy landscape shifts to 



nonzero polarization or magnetization. 

Quantitatively, to study the FMPT in dilute Fermi gases at 
zero temperature, we should calculate the energy density & 
as a function of the spin polarization or magnetization x = 
(«f - «|)/(«-j- + ni) at gi yen dimensionless gas parameter k?a 
which represents the interaction strength H ■ Here, fcp is the 
Fermi momentum related to the total density n = w-r+nj, by n — 
fcp/(37r 2 ) and a > is the s-wave scattering length. Generally, 
the energy density can be expressed as 



8(x) = -nE F f(x), 



(1) 



where E-p = k 2 /(2M) is the Fermi energy with M being the 
fermion mass. The dimensionless function f(x), which de- 
pends on the gas parameter k F a, represents the energy land- 
scape with respect to the magnetization x. 

For the order of the FMPT, Belitz, Kirkpatrick, and Vojta 
(BKV) [3j] have argued that the phase transition in clean itin- 
erant ferromagnets is generically of first order at low tempera- 
tures, due to the correlation effects or the coupling of the order 
parameter to gapless modes that lead to a nonanalytic term in 
the free energy. The general form of the Ginzburg-Landau 
free energy for clean itinerant ferromagnets takes the form 



/glO) = ax 2 + ux 4 ln|jc| + fix 4 + 0(x 6 ), 



(2) 
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where we can keep /3 > 0. If the coefficient v is positive, the 
phase transition is always of first order. On the other hand, for 
negative v, one always has a second order phase transition. 
The BKV argument is based on the assumption v > moti- 
vated by perturbation theory ll37ll . This is true for many solid 
state systems where the FMPT occurs at weak coupling. How- 
ever, for dilute Fermi gases where the critical gas parameter is 



2 



expected to be of order (9(1), the assumption of a positive v is 
not reliable. 

In this paper, we will study the non-perturbative effects on 
the FMPT by summing a set of particle-particle ladder dia- 
grams to all orders in the gas parameter, motivated by the 
large-dimension expansion proposed by Steele HH. We con- 
sider a universal repulsive Fermi gas where the effective range 
effect can be neglected, corresponding to a two-component 
"upper branch" Fermi gas with a positive s-wave scattering 
length. The prediction may also apply to the hard sphere gas 
since the effective range corrections are subleading order con- 
tributions in the large-dimension expansion. Our main conclu- 
sions for the order and the critical gas parameter of the FMPT 
can be summarized as follows: 

(1) Order of phase transition. We predict a second order phase 
transition, in contrast to the BKV argument. This suggests that 
the FMPT in dilute Fermi gas may correspond to the case of 
negative v. 

(2) Critical gas parameter. We predict a critical gas parameter 
(kpd) c = 0.858 where the spin susceptibility \ diverges. The 
critical gas parameter and the spin susceptibility we obtained 
are in good agreement with the Quantum Monte Carlo results 

a. 

The paper is organized as follows. In Section II we briefly 
review the perturbative predictions for FMPT in dilute Fermi 
gases. In Section III we introduce the effective field theory ap- 
proach to the two-body scattering problem and show how we 
can recover the scattering amplitude by ladder resummation. 
We study the non-pertuabtive effects on FMPT in the theory 
of ladder resummation in Section IV and investigate the role 
of hole-hole ladders in Section V. We summarize in Section 
VI. 



II. PERTURBATIVE PREDICTIONS 

In the perturbation theory, the gas parameter kpa is treated 
as a small number. Up to the order 0((kpa) 2 ), the expression 
for f(x) is universal, i.e., independent of the details of the 
short range interaction. We have 



1 



f(x) = -(^ + 7$ + 



10k F a 
~9n~ 



3„3 



(ha) 2 
21n 2 



f^T-^). ( 3 ) 



where r/i = (1 + x) 1 ^ 3 and rji = (1 — x) 1 ^ 3 . The Oth-order 
term corresponds to the kinetic energy, and the lst-order term 
coincides with the Hartree-Fock mean-field theory The 
coefficient ((jli, ?7j.) in the 2nd-order term was first evaluated 
by Kanno IB9I1 . Its explicit form is 

% = 2277^(77 T + ?7 i )- 477{lii 477[ln- 



'7r 
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(4) 



Setting x — 0, we recover the well-known equation of state for 
dilute Fermi gases, 



1 10 7 

1 H kva 

9n 



4(11 -21n2) 
2ht 2 



(k F aY 



(5) 



which was first obtained by Huang, Yang, and Lee 14011 and 
recovered by Hammer and Furnstahl [PlU in recent years using 
effective field theory. 

In the lst-order PTh, the FMPT is of second order and 
occurs at k^a = n/2 J^]. However, taking into account 
the 2nd-order corrections, one finds a first order FMPT at 
kpa = 1.054 oi m. This can be understood by noticing 
the non-analytical term <x x 4 ln|x| with positive coefficient in 
the small-x expansion of the coefficient £(?7f , 77 j,). The small-x 
expansion of the energy density (3) takes the form 

f( x ) = /(0) + ax 2 + vx 4 \n\x\ + fix 4 + 0(x 6 ), (6) 

which is consistent with the BKV argument i37tl that the cor- 
relation effects or the coupling of the order parameter to gap- 
less modes generally leads to non-analytical terms in the free 
energy. The coefficient u can be evaluated as 



40(k F a) 2 
243tt 2 



(7) 



Therefore, up to the order 0((k F a) 2 ), the Fermi gas problem 
corresponds to the case v > which is assumed in the BKV 
argument. 

In general, we expect that the critical parameter is of order 
0(1). Therefore, the perturbative predictions for the FMPT 
are probably unreliable. There naturally arises a serious prob- 
lem: Does the dilute Fermi gas problem really correspond to 
the case u > if the non-perturbative effects at k?a ~ 0(1) are 
taken into account? For the two-body problem in the vacuum, 
it is well known that an infinite set of bubble diagrams with the 
leading-order contact interaction must be resummed in order 
to reproduce the correct scattering amplitude if the two-body 
scattering length is large ll42ll . Therefore it is natural to extend 
the idea of ladder resummation to finite density so that the pre- 
dicted equation of state works well even at kpa ~ 0(1). We 
can also compare the non-perturbative predictions with the re- 
sults from recent Quantum Monte Carlo (QMC) simulations 



m. LADDER RESUMMATION FOR TWO-BODY 
SCATTERING 

Before we establish a non-perturbative description for the 
FMPT in dilute Fermi gases, it is instructive to start with the 
low energy effective field theory in vacuum ll4lll42ll and to see 
how the two-body scattering amplitude is reproduced from the 
ladder resummation method. 

For non-relativistic two-body scattering in the s-wave chan- 
nel associated with a short range interaction, the scattering 
amplitude J[(k) is related to the s-wave scattering phase shift 
6 by 



m.(k) 



An 



1 



M k cot 6 - ik 



(8) 
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where k = |k| is the scattering momentum in the center-of- 
mass frame. If there exist bound states for attractive interac- 
tions, the scattering amplitude should exhibit some imaginary 
poles, k = i -\J-ME\„ on the complex k plane with E\, < 
being the binding energy. In general, the short range interac- 
tion is characterized by a momentum scale A. Therefore, for 
low energy scattering, i.e., k <k A, the quantity kcotd can be 
expanded as a Taylor series in k 2 /A 2 . In quantum scattering 
theory this is called the effective range expansion. 



kcotd 



1 1 ,2 

a 2 



(9) 



where a is the scattering length, and ro is the effective range. 
For a natural system [42], we have \a\ ~ 1/A and \r„\ ~ 1/A. 
An example commonly studied is a hard-sphere gas with ra- 
dius R, in which case a - R and ro = 2/?/3. For cold atomic 
gases the inter-atomic interaction can be tuned by means of 
the Feshbach resonance, we can have \a\ » \r„\ ~ 1/A. 

According to the effective range expansion, one can con- 
struct the low energy effective field theory ll4ll l42ll describing 
scattering at momenta k <K A. Since we assume k <K A, all in- 
teractions in the effective Lagrangian are contact interactions. 
The low energy effective Lagrangian contains infinite contact 
interaction terms and is given by fiHl 



-Ceff = W 

c 2 

+ — 

16 



id, 



2M 

(i/^) 1 " (^Vfy) + h.c 



(10) 



where Co and C2 are dimensionful coupling constants and 
<->—><— 

V = V - V is a Galilei invariant derivative, and . . . denotes 
interactions with more derivatives (~ V 2 ", n > 2) which gen- 
erally have coupling constants Cn n . The coupling constants 
C2 n (n = 0, 1, 2, ...) should be determined by reproducing the 
scattering amplitude tft{k). 

In practice, we can reproduce the scattering amplitude Jttk) 
order by order in a Taylor expansion in k/A. For small scat- 
tering length (\a\ ~ 1/A and \ak\ <K 1), we can expand the 
scattering amplitude as 



4na 
m) = ~~M 



1 - iak 



(ar 2 \ 2 



(id 



However, for large scattering length (\a\ » 1/A), Kaplan et 
al. showed that one needs to expand 3\kk) in powers of k/A 
while retaining ak to all orders II42II : 



An 1 
M 1 /a + ik 



r /2 
l/a+ ik' 



■k 2 



(12) 



This means, if the scattering length is large, the loop diagrams 
with the leading order interaction Co have to be resummed. 

According to the free fermion propagator @o(po, p) = 
l/(po _ w p + ie) with the free dispersion oj p = p 2 /(2M), the 
one-loop bubble diagram Bq(Pq, P) (FigQJ) is given by 



BoCPo,P) 



. f d 4 g 

J (2tt) 4 ^ + 



9o 



(P/2+q) 2 • Po 
2M T ,fc 2 



qo 



(P/2-q) 2 

2M 



+ ie 



r d\ 1 



(13) 



Here Po and P are the total energy and momentum of the pairs 
in the bubble diagram (see FigQla)). Let p! and P2 be the 
momenta of the scattering fermions, we have P = pi + P2 
and k = (pi - pi)/2. Further, if the on-shell condition Po = 
(Pi + P?)/( 2M ) = P 2 /(4M) + k 2 /M is imposed, we find that 
Bo depends only on the relative momentum k, corresponding 
to the translational invariance. 

The integral over q in Eq. ( TT31 is linear divergent and there- 
fore needs to be regularized. A natural regularization scheme 
is to use a momentum cutoff equal to A 14311 . In this paper, we 
employ the dimensional regularization scheme. To this end, 
we change the space-time dimension from 4 to D and multi- 
ply the integral by a factor (/j./2) D . Here fi is an arbitrary mass 
scale introduced to allow the couplings Co,, multiplying oper- 
ators containing V 2 " to have same dimensions for any D. In 
general, the integral Bo in D dimension can be evaluated as 



B„(p ,P) - -r^— j^zr^ 

/ p2 \(P-3)/2 

x Ml-MP + — -ie\ . (14) 

For small scattering length, it is convenient to use the minimal 
subtraction (MS) scheme which subtracts any l/(D - 4) pole 
before taking the D — > 4 limit. However, for large scatter- 
ing length, it is more convenient to use the power divergence 
subtraction (PDS) scheme. The PDS scheme involves sub- 
tracting from the dimensionally regularized loop integrals not 
only the 1 /(D - 4) poles corresponding to log divergences, as 
in MS, but also poles in lower dimensions which correspond 
to power law divergences at D = 4. The integral Bo has a pole 
in D = 3 dimensions. It can be removed by adding a counter 



4 



term 6B = Mp/[An(3-D)] to B 14: 
integral in D — 4 dimensions is 



Finally, the subtracted 



B (P ,P) = ~ 
An 



p2 

-MP + — -ie 



(15) 



Note that the MS scheme corresponds to the p = case. 

The dependence of C2„(a<) on /" is determined by the re- 
quirement that the scattering amplitude is independent of the 
arbitrary mass scale p. To this end, we impose the on-shell 
condition, P = P 2 /(4M) + k 2 /M. Then the one-loop bubble 
diagram becomes Bo(k) = -M(p + ik)/(An). Summing the 
bubble diagrams with Co vertices, we obtain ll42tl 



Coin) 



Ci(p)k 2 



1 - C (ti)Bo(k) [l - C (/i)Bo(^)] 
Comparing this result with the expansion (fT2] i, we obtain 
An 1 



(16) 



C (ju) 



C2(M) = £ 



M -p+ 1/a' 
1 



M \-fi + 1/a 



2 ' 



(17) 



It was shown that these results fulfill the renormalization 
group equations ll42ll . We note that the mass scale p is sim- 
ilar to the cutoff A. In the cutoff scheme, we have Co(A) = 
(An/M)(-2A/n + 1/a) 0. 

In the following, we mainly consider a short range poten- 
tial with a positive scattering length a and negligible effective 
range ry <*c a. In this case, we are able to obtain a univer- 
sal result for f(x) which is independent of the details of the 
interaction. In this case, the pair propagator Sq(Pq, P) in the 
vacuum is given by 



S {P ,P) = 



Co(yu) 



l-C p (?o,P) 
An 1 



M 



1/a- yj-MPo + f 



(18) 



For positive scattering length, the pair propagator has a pole 
given by Pq = -l/(Ma 2 ) + P 2 /{AM). This pole corresponds 
to a bound state with binding energy E h = -I /(Ma 2 ) and 
effective mass 2M. Therefore, if the effective range is neg- 
ligible, the underlying potential must be attractive and the 
ground state is a bound molecule of size a. However, for 
two-body scattering state with positive center-of-mass energy 
E = k 2 /M > 0, the effective force is repulsive. This is the so- 
called "upper branch" which is well defined in the two-body 
picture. For the many -body problem, a metastable "repulsive" 
Fermi gas can be realized if all fermions are forced on the up- 
per branch of a Feshbach resonance with a positive s-wave 
scattering length lHH. 



IV. LADDER RESUMMATION AT FINITE DENSITY 

We now turn to the many-body problem of the repulsive 
Fermi gases. The main purpose of this paper is to present 



(a) 



(b) 




Pi 



B 



FIG. 1: (a) The elementary particle-particle bubble B(p^,p{) with 
external momenta p\ and p2 for the two spin components. The solid 
line with arrow corresponds to the particle term of the propagator 
H9\ . The dashed line represents the interaction vertex Co- (b) A 
typical particle-particle ladder diagram contributing to the interaction 
energy. 



a non-perturbative calculation for function f(x) from which 
we can give a better prediction for the FMPT In general, we 
expect that the new non-perturbative result for f(x) satisfies 
the following two criteria: (i) The function f(x) recovers the 
perturbative result Eq. (|3) at weak coupling kpa — > 0; (ii) 
Since we consider a short range potential with ro <K a, the 
physical result should be universal, i.e., f(x) depends only on 
the gas parameter kpa and does not depend on other param- 
eters such as the renormalization scale fi. The criterion (ii) 
is hard to be fulfilled since the loop corrections in quantum 
field theory generally brings the renormalization scale depen- 
dence and weakens the prediction power due to the arbitrari- 
ness in the choice of the renormalization scale. However, in 
the following we will show that the result from the particle- 
particle ladder resummation, corresponding to the leading or- 
der of the large-dimension expansion, is independent of the 
arbitrary mass scale pi. 



A. Non-Perturbative Energy Density 

We first construct the non-perturbative version of the energy 
density S(x) using the vertex Co(^) determined in Section III 
and the free propagators for the two spin components at finite 
density [45] 

Q(|p|-*g) gCg-lgj) 

Saipo, P) = V + ~ -■ (19) 

po - u> p + ie po - a>p- ie 

Here cr =f , J,, — kprj-^ i are the Fermi momenta of the two 
spin components, and 0(z) is the Heaviside step function. For 
each spin component, the propagator (fT9b describes two types 
of excitations, particles with momentum |p| > k^ and holes 
with |p| < k%. 

The dilute imperfect Fermi gases are best described by re- 
summing the multiple interactions in terms of the scatter- 
ing amplitude. The Galitskii integral equations [45] for the 
effective two-particle interaction or scattering amplitude in 
the medium are given by the ladder resummation ll46ll . On 
the other hand, for large gas parameter k^a, one may look 



5 



for other expansion parameters instead of k-pa itself. Steele 
ll38tl and Schafer et.al. i47ll have suggested a new expan- 
sion method using l/D as the expansion parameter, where 
J) = 2 D I 2 with D being the space-time dimension. Most 
importantly, they have shown that the contribution of the 
particle-particle (pp) ladder resummation fiff , is the leading- 
order contribution of the 1 ID expansion 11381 14711 . i.e., 

& = &^ n +&^ + 0(l/D). (20) 

All other contributions like hole-hole (hh) ladder sum and ef- 
fective range corrections are suppressed by a factor 1 /D. Ac- 



cording to the above arguments, we expect that the most im- 
portant non-perturbative contributions come from the leading 
order of the 1/D expansion. The interaction energy density 
fijjj' to this order is given by all the particle-particle scattering 
terms (i.e., the npp-lhh bubbles for all n — 0, 1, 2, • ■ ■ ). 

To evaluate the interaction energy density, we first calculate 
the elementary in-medium particle-particle bubble B(Po,P) 
shown in Fig. [TJa). The fermion lines in the bubble dia- 
gram correspond to the particle terms of the free propagator 
(fT~9b . According to the finite density Feynmann rules B41H . it 
is given by 



B(P ,P) 



. r d*g W 
J (2tt) 4 ^ + 



0(|P/2 + q| 



0(|P/2-q|-4) 



qo 



(P/2+q) 2 • Po 
2M T ,fc 2 



qo 



(P/2-q) 2 
2M 



d\ ®(|P/2 + q| - 4)0(|P/2 - q| - k l F ) 



{2nf 



- 21 + ie 

AM M 



(21) 



For vanishing densities, = 0, the in-medium particle-particle bubble recovers the vacuum result Bo- If the on-shell condition 
is imposed, the in-medium particle-particle bubble B depends on not only the relative momentum p but also the total momentum 
P. This is due to the loss of translational invariance in the presence of Fermi sea. 
We can separate B into a vacuum part and a medium part using the identity 

0(|P/2 + q| - 4)©(|P/2 - q| - 4) = 1 - ®(k\ ~ |P/2 + q|) - 0(4 - |P/2 - q|) + 0(4 - |P/2 + q|)0(4 - |P/2 - q|). (22) 

The vacuum part (corresponding to 1) is identical to Bo defined in the last section and is linearly divergent. The medium part 
is convergent. For the vacuum part, it is natural to use the dimensional regularization with PDS scheme introduced in the last 
section. 

Then the zipp-lhh bubble (see Fig. [TJb) for a typical example) at given n reads 



4 k 



0(4 
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+ k ° 2M e T 



> (P/2-k) 2 ■ 
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" ~ ° J (2n)* J (2n 
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where e"' p ° with r/ — > + is a convergence factor 114 ill . The 
integration over Po picks up the pole or imposes the on-shell 
condition P = P 2 /(4M) + k 2 /M, i.e., 



dP j nPo [B(P ,V)r 
2ni 6 Po - — 



M Lt 



[fi(P,k)]'\ 



where the on-shell version of B is given by 
B(P, k) = 

d\ 0(|P/2 + q| - 4)©(|P/2 - q| - 4) 



M 



J (2nf 



k 2 - q 2 + ie 



The total interaction energy density 6[? is given by 



n=0 



(24) 



(25) 



(26) 



Completing the summation of this geometric series, we obtain 
the interaction energy density at the leading order of the 1/2) 
expansion, 



o(0) = c f d Pi f d & 
int °J (2tt) 3 J {2nf 



d 3 p 2 0(4 -|pi 1)0(4 -|p 2 |) 



l-Cofi(P.k) 



, (27) 



where pi ; 2 = P/2 ± k as defined in Section EH. The imaginary 
part of B can be evaluated as 



EmB(P,k) = -^0(|pil -4)0(|p 2 | -4)- 



(28) 



This quantity is nonzero only when the momenta pi and p2 
are both above the Fermi surfaces. However, the final inte- 
gration over p! and P2 in the interacting energy density £i nt is 
associated with a phase-space factor 0(4 ~ IPi 1)0(4 ~ Ip2l)- 
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Therefore, the interaction energy density is real and physical, 
as we expected. 




FIG. 2: The pole energy E (divided by the binding energy |£t,l = 
1 1 {Ma 1 ) in the vacuum) at zero pair momentum P = as a function 
of the gas parameter kpa. The pole energy turns out to be positive for 
kpa > n/2. 




FIG. 3: (Color-online) The pole energy E (divided by 2Ef) as a func- 
tion of the pair momentum P = |P| (divided by kp) for various values 
of the gas parameter k r a. The dashed line corresponds to the disper- 
sion E(P) = P 2 I(AM). 



B. In-Medium Two-Body Problem 

Since we adopt a zero range potential, a bound state with 
binding energy E\> = -1 /(Ma 2 ) always exists for positive 
scattering length a > 0. A key problem here is that how 
we can describe a metastable repulsive Fermi gas where all 
fermions are forced on the scattering states. Actually, the so- 
called "upper branch" has clear meaning only in the two-body 
picture, and so far it is not clear to what extent this two-body 
picture of a "repulsive" Fermi gas will persist. Recent study 



on three attractive fermions shows that there are many non- 
trivial avoided crossings between the two branches close to 
the resonance (a — > oo), making it difficult to unambiguously 
identify a repulsive Fermi system lfl7ll . 

To realize a metastable repulsive Fermi gas we have to ex- 
clude the molecule bound states of two atoms with unlike 
spins and enforce all atoms to the scattering states 114 811 . One 
possible prescription is to subtract the contribution from the 
bound-state poles within the Nozieres-Schmitt-Rink (NSR) 
theory ll48tl . However, as designed, NSR theory works well 
only at temperature higher than the critical temperature of 
superfluidity. The particle-particle resummation theory we 
present complements the NSR theory and can be regarded as 
the zero-temperature analogue of the NSR theory. 

The key point in this problem is to consider the medium 
effects on the bound state properties. To this end, we first con- 
struct the pair propagator S(Pq, P) in the presence of Fermi 
seas. With the in-medium elementary particle-particle bubble 
B(Pq, P), the in-medium pair propagator S(Po, P) is given by 
the ladder resummation, 



S(P ,V) = 



CoOu) 



I - C (m)B(Po,P) 



(29) 



With this pair propagator, the interaction energy density £ 
can be expressed as 



(0) 



? (0) 



r d s ~p r 



IPll)®(4-|P2l) 



J 2^ 



S(P ,P) 



"0 — 777 — 77 



(30) 



In general, the in-medium pair propagator S(Pq, P) has a 
real pole Pq = E(P) corresponding to the in-medium bound 
state. However, we now show that such pole does not con- 
tribute to the energy density for the regime of the gas param- 
eter k F a we are interested in. As shown in the last subsection, 
in the calculation of the interaction energy density <5j nt , the 
on-shell condition Pq = P 2 /(4M) + k 2 /M is imposed and the 
integrations over the momenta P and k are performed accord- 
ing to the finite density Feynmann rules. Therefore, if the 
energy dispersion of the pole E(P) satisfies the condition 



p2 

EfP) < 

V ; 4M 



(31) 



for arbitrary P, its contribution to the energy density is natu- 
rally excluded. 

Since the main purpose of this paper is to study the 
FMPT which corresponds to an instability toward a small 
polarization x, we can set x = here. For convenience, 
we define two dimensionless quantities s = \P\/(2kp) and 
z = y/MP - P 2 /4 + ie/kp = ^P /(2E F ) - s 2 + ie. The in- 
medium pair propagator can be evaluated as 



S(P ,P) 



An 



1 



M l/a-(k F /jT)W(s,z) 



(32) 



where W(s, z) is given by 
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W(s,z) = 



1+s-z 1 
1 + s + z In h 



1 + s + z 



s 2 -z\ (l + s) 2 -z 
In 



2s 



(l-z) 2 -s 2 l-/-z 2 , 
2 + z In — — — + = In 



1 - s 2 - z 2 
(l+s) 2 -z 2 



s 2 -(l +z) 2 



2s 



(l-s) 2 -z 2 



©(1 - s) 

0(5-1). 



(33) 



The same result was also obtained in a recent paper ll49ll . 

For zero pair momentum P = 0, the condition £(P) < 
P 2 /(4M) implies £(0) < 0. We thus focus on the regime of 
the gas parameter where the pole £(0) is negative. In this case, 
£(0) is determined by a simple equation, 

7T, — = 1 + a/ -TE— arctan V ~TE~- (34) 
2k F a V 2£ F V 2£ F 

This equation has negative solution only for < kpa < n/2 
where the solution represents the binding energy of a in- 
medium bound state. The numerical result for £(0) is shown 
in Fig. [2] In the low density limit kp — ♦ 0, £(0) recovers 
the vacuum result £ vac (0) = £\, = -I /(Ma 2 ). However, at 
finite kp, the medium shields the bound state and reduces the 
binding energy, i.e., |£(0)| < |£ V ac(0)l' For kpa > n/2 and 
kpa < 0, Eq. (f34-b has a positive solution which corresponds 
to the positive energy pole of the in-medium pair propagator. 
Such a pole is associated with Cooper pairs and its appear- 
ance represents the BCS instability. This positive energy pole 
does not lead to singularities in the energy density integration, 
as can be seen in next subsection, and does not need special 
treatment. 

Note that £(0) < is not a sufficient condition for £(P) < 



P 2 /(4M). We thus have to check the energy dispersion £(P) 
carefully. The numerical results for some values of the gas pa- 
rameter kpa are shown in Fig. [3] For kpa < 1 .34, the condition 
£(P) < P 2 /(4M) is fulfilled for all values of P. However, for 
1.34 < kpa < n/2, there exists a regime P\ < |P| < £2 where 
£(P) > P 2 /(4M). 

In conclusion, the condition £(P) < P 2 /(4M) is fulfilled for 
kpa < 1.34. Therefore, in the parameter regime kpa < 1 inves- 
tigated in the following, the contribution from the bound state 
can be naturally excluded in the ladder resummation scheme. 



C. Evaluating the Energy Density 

Now we evaluate the explicit form of the energy density 
£(x) and the dimensionless function f(x). First, the elemen- 
tary particle-particle bubble B(P, k) can be decomposed into 
four parts 

B(P, k) = fl (P, k) + fl T (P, k) + B X (P, k) + %(P, k), (35) 

where £0 is the vacuum part discussed in Section III, and the 
other parts are given by 



d\ 0(4 - |P/2 - q|) 
ie 



C d'q ®(k' - P/2 + q) r d'q &(k'- P 2 

J (2n) 3 k 2 - q 2 + ie J (2n) 3 k 2 - q 2 + 

/xj ■ \ I# r d\ ml - IP/2 + q|)0(4 - |P/2 - q|) 

B n (P,k) = M — — 2 — . ■ 

J (2ny k 2 -q 2 + ;e 



(36) 



For convenience, we define another dimensionless quantity t = |p|/^F together with s defined in the last subsection. Since the 
imaginary part of B does not contribute to the interaction energy, we need only to evaluate the real part of B. We have 



Ma Mkp Mkp Mk F 

ReB (s, t) = RsBf(s, t) = j^R^ 0, ReB A (j, = f), ReB n (s, t) = ^f%(^ 0- 



where Ro-(s, t) (cr =f , J,) reads 

R a (s,t) = 
and the function Rfi(s, t) is 

Rn{s,t) = \ 



■ (s + ty 



4s 



In 



TJa- + S + t 



S - t 



ril-(s-tf 
4s 



In 



TJtr + S-t 



s + t 



+ llo-, 



-Q(x)R l (s, t) - ®(-x)Rj(s, t) , < s < i|77 T - 77J 

K 1 (s,t) + K L (s,t) , i|77 T -77J < s < i|j/ T +774,1 
, elsewhere. 



(37) 



(38) 



(39) 
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Here A"o-(s, t) is defined as 



As 



-In 



for - ^) 2 - f 2 



r 2 _ ,2 _ ,2 



+ 2 ln 



77o- - 5 + t 



S-TJ 



- S-t 



(40) 



where r = (77^ + ?7 2 )/2. 

Finally, the elementary particle-particle bubble reads 

Mu Mk F 

B(s,t) = + -^-R w {s,f), (41) 

where the function R pp (s, t) is defined as 

R pp (s, t) = R t (j, f) + f) + « tt (j, 0, (42) 

Substituting the result of B(s, t) into the expression of we 
observe that the energy density is independent of the renor- 
malization mass scale fi. Converting the integration variables 
Pi and P2 to P and k, we find that the function f(x) can be 
expressed as 

1 so f °° r°° 

fix) = ^ + — J * 2 ds J Q tdtI(s,t)F(s,t), (43) 
where F(s, t) is given by 
F(s,t) = 



k F a 



1 - j-k F aR pp (s, t) 



(44) 



The function I(s, t) appears due to integration over the angle 
between P and k. Its explicit form is 



I(s,t) = 



tf-(s + t) 1 



As 



-®(s + t-7]<t) + (t7 T -> 774.) + t 
f 2 )0(/7 T -| S -f|)0(^-k-f|). (45) 



As we mentioned in the beginning of this section, it is im- 
portant to check whether the present result for f(x) is con- 
sistent with the perturbative expression (|3) for weak coupling 
kpa <k 1. To this end, we expand the function F(s, k) as 

F(s, t) = k F a + -(k F a) 2 R pp (s, t) + 0((k F a) 3 ). (46) 

71 

Using the expressions for I(s, f) and R pp (s, t), we can show 
that 

so r°° , r°° 10 , , 

— s z ds tdtl(s, t) = — 77^77; (47) 
n Jo Jo 9?r 1 



and 



so p 2 r 



tdtl(s, t)R pp (s, t) = 



2ln 2 ' 



(48) 



Therefore, our non-perturbative expression (27) exactly re- 
cover the perturbative result (3) at weak coupling. This con- 
vinces us that the present theoretical approach is suitable to 
study the universal upper branch Fermi gas with a positive 
scattering length. In addition, we can compare the results from 
our theory and the second-order perturbation on the same foot- 
ing and study the non-perturbative effects on the FMPT 



D. Results and Discussion 

(A)Energy density and compressibility. We first study the 
equation of state for the unpolarized case x = 0. The gas 
parameter dependence of the energy density & in the regime 
< k F a < 1 are shown in Fig. |4] We find that the result from 
the ladder resummation is consistent with the perturbative re- 
sult (3) for small gas parameters k F a < 0.4. However, signifi- 
cant deviations are found for k F a > 0.4, consistent with recent 
Quantum Monte Carlo simulations fHIEl]. For the Quantum 
Monte Carlo simulations of the attractive interactions with a 
negligible effective range (corresponding to UB and UB2 in 
Fig. HJ, the exclusion of molecular bound states is imple- 
mented by choosing a two-body Jastrow factor JH |^] to be 
the scattering solution of the attractive potential correspond- 
ing to positive energy, which, by construction, is orthogonal 
to the bound molecules. Therefore, the accuracy of the Quan- 
tum Mote Carlo data depends on the choice of the Jastrow 
factor. Actually, exact orthogonality of the many-body vari- 
ational wave function to the superfluid ground state (molecu- 
lar condensation) can not be achieved in the Quantum Monte 
Carlo simulations We note that our theoretical curve 

agrees better with the UB2 data than with the UB data. The 
reason could be that the UB2 data from [6J are obtained with 
a Jastrow factor which imposes a better orthogonality to the 
superfluid ground state. 

An important issue is whether the system is mechanically 
stable. The mechanical stability of the system requires a pos- 
itive compressibility a:, which is defined as 



1 2 d 2 S 
k on- 



(49) 



For the present ladder resummation theory, the explicit form 
of k can be evaluated as 

— = 1 + I s 2 ds tdtl(s, t)G(s, t), (50) 

« n Jo Jo 

where kq = 3/(2nE F ) is the compressibility for non- 
interacting Fermi gases, and the function G(s, t) is given by 



G(s, t) = F(s, f) + 



5R pp (s,t)^ 2/ _ flpp(M)^ 



9n 



-F z (s,t) + 



9jt 2 



F 3 (s,t). (51) 



The compressibility k as a function of the gas parameter k F a 
is shown in Fig. [5] Comparing to the result from the second- 
order perturbation theory, 

ko ,2, 8(ll-21n2), f . 2 

— = 1 + -k F a + — (k F a)\ (52) 

k n 157r z 

good agreement is found for small gas parameters, as we ex- 
pected. In the regime < k F a < 1 we are interested in, we 
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FIG. 4: (Color-online) The energy density £ (divided by its value 
Go = |;j£*f f° r non-interacting Fermi gases) as a function of the gas 
parameter kpa (0 < k^a < 1) for the unpolarized case x = 0. The 
solid line is the result calculated from our particle-particle ladder 
resummation theory (RTh). The dashed line is result of the 2nd- 
order perturbation theory (PTh). The dash-dotted horizontal line 
corresponds to the energy of the fully polarized state (x =1), i.e., 
f(l) = 2 2,i . The blue squares are the Quantum Monte Carlo (QMC) 
data for the hard sphere (HS) potential @|, the red circles for the 
upper branch (UB) of a square well potential j5|], and the green di- 
amonds for the upper branch (UB2) of an attractive short range po- 
tential For UB and UB2 cases, the effective range r is much 
smaller than the s-wave scattering length a (HQ]. 



"5 




FIG. 5: The compressibility k (divided by its value ko = 3/(2nEp) 
for non-interacting Fermi gases) as a function of the gas parameter 
k^a (0 < kpa < 1) for the unpolarized case x = 0. The solid line 
is the result calculated from our particle-particle ladder resummation 
theory. The dashed line is result of the 2nd-order perturbation theory. 



find that the compressibility k is positive, indicating that the 
system is mechanically stable. 

(B)Spin susceptibility. Next we study the response of the 



energy density to an infinitesimal polarization x. This re- 
sponse is referred to as the spin (or magnetic) susceptibility. 
The spin susceptibility x can be defined as 



I - 1 d 2 £ 

X n 2 dx 2 



3E F d 2 f(x) 



x=0 



5n dx 2 



(53) 



A=0 



In the present ladder resummation theory, an explicit form of 
X is hard to obtain. In practice, we expand the function f(x) 
near x = as f(x) = f(0) + ax 2 + ■ ■ ■ . The coefficient a is 
related to the spin susceptibility by 



Xo 
X 



(54) 



where xo = 3n/(2£p) is the spin susceptibility of non- 
interacting Fermi gases. Therefore, a diverging spin suscepti- 
bility generally indicates a FMPT, as long as the transition is 
of second order. 

In the second-order perturbation theory, an analytical result 
for^- can be achieved, 



x 



l 



> 16(2 + ln2) 2 

KpU \K-F U ) I 

r 15n~ 



(55) 



which indicates a diverging spin susceptibility at kpa = 1.058. 
However, this differs from the critical gas parameter {kpa) c = 
1.054, because the phase transition is of first order in the 
second-order perturbation theory due to the appearance of the 
non-analytical term vx A In \x\ with u > 0. 

Our result for the spin susceptibility ^ as a function of the 
gas parameter kpa is shown in Fig. [6] and compared with the 
perturbative result. We find that the spin susceptibility pre- 
dicted by the ladder resummation deviates significantly from 
the second-order perturbative result for kpa > 0.4. Further, 
the spin susceptibility diverges at kpa = 0.858, in contrast 
to the value 1.058 from the second-order perturbation theory. 
The data from the Quantum Monte Carlo simulations JH are 
also shown in Fig. [6] as a comparison. Our theoretical re- 
sult is in good agreement with the data for the upper branch 
of the square well potential where the effective range ro is 
tuned to be much smaller than the scattering length H. The 
gas parameter kpa = 0.86 where x diverges is very close to 
our prediction kpa = 0.858. For the purely repulsive poten- 
tial, i.e., the hard sphere potential, the effective range effect 
can not be neglected in prior. However, we find that our re- 
sult still has nice agreement with the data for the hard sphere 
case. The gas parameter kpa = 0.82 where x diverges is also 
close to our prediction k F a = 0.858. Actually, the differ- 
ence between the upper branch and the hard sphere cases, i.e., 
0.86 - 0.82 = 0.04, is very small compared with the critical 
gas parameters. This indicates that the contribution from the 
effective range effect is relatively small even for k-pa ~ 0(1), 
if the Quantum Monte Carlo results are reliable. This can be 
understood from the large-dimension expansion ||3~8"[ l47ll in- 
troduced in the beginning of this section: The particle-particle 
ladder sum is the leading-order contribution in the 1 /£) expan- 
sion, and all other contributions including the effective range 
corrections are suppressed by a factor 1 I'D. 
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FIG. 6: (Color-online) The dimensionless inverse spin susceptibility 
Xolx as a function of the gas parameter k f a. The blue squares and red 
circles are the QMC data H 

for the HS and UB cases, respectively. 
The solid line is the result calculated from the particle-particle ladder 
resummation. The dashes line is the 2nd-order perturbative result. 



turbation theory 14111. However, some definite conclusions can 
be drawn from our numerical results: (1) Higher-order terms 
in the gas parameter can also generate non-analytical terms 
of the form x 4 ln|x| and may generate other important non- 
analytical terms which are not known due to the mathemati- 
cal limitation; (2) The coefficients of the non-analytical terms 
generated by the higher-order contributions are certainly not 
always positive, and they are generally proportional to (kpa)" 
for the n-th-order contributions. Since the phase transition 
occurs at a gas parameter kpa ~ 0(1), the non-perturbative 
effects from the sum of the higher order contributions are very 
important. As we have shown numerically, their effects are 
not only reducing the critical value of the gas parameter but 
also changing the order of the phase transition. 



(C)Ferromagnetic transition. While a diverging spin sus- 
ceptibility indicates a ferromagnetic phase transition, the or- 
der of the ferromagnetic phase transition and the critical gas 
parameter (kpa) c should be obtained by studying carefully 
the shape of the energy landscape, i.e., the full jt-dependence 
of the function f(x). To very high numerical accuracy, we 
haven't found any maximum at x + in the energy landscape. 
Instead, we find a second orderphase transition at kpa = 0.858 
where the function f(x) starts to develop a minimum at x + 0, 
consistent with the gas parameter where the spin susceptibil- 
ity diverges. This is in contrast to the second-order pertur- 
bation theory which predicts a first order phase transition at 
kpa = 1 .054 ifTill where the spin polarization x jumps from 
zero to x c = 0.573. A second order FMPT for a zero range 
potential model was also obtained by Heiselberg IU8I1 recently 
using a completely different many-body method. 

It seems that our result of a second or der p hase transition 
is in contradiction to the BKV argument 0711 . However, the 
BKV argument is based on the assumption that v > 0. Ac- 
tually, we have fitted the energy density of the form f(x) = 
f(Q) + ax 2 + vx 4 ln\x\ + fix 4 for small x. For small gas param- 
eter kpa < 0.3, the coefficient u agrees well with the pertur- 
bative result v = 40(A:Ffl) 2 /(2437r 2 ). However, for larger kpa 
(especially around the critical gas parameter), it turns out to 
be negative due to the non-perturbative effects. This indicates 
that the FMPT in the systems of dilute repulsive Fermi gases 
corresponds to the case v < 0, and is a counterexample to the 
BKV argument where the assumption v > is adopted. 

Since an analytical expression for the function f(x) as well 
as the coefficient v cannot be achieved in the present ladder re- 
summation theory, we cannot understand analytically how the 
non-partuebative effects modify the order of the phase tran- 
sition. In fact, analytical results cannot be obtained from the 
order 0((kpa)^) even for the unpolarized case x — in the per- 



FIG. 7: The elementary bubbles organized in the number of the MI. 
The solid line with a cut represents the MI part of the propagator, and 
the pure solid line corresponds to the vacuum part. 



V. INCLUSION OF HOLE-HOLE LADDERS 

In this section we check whether our conclusion that the 
FMPT is of second order is changed by other contributions. 
We consider the contributions from the hole-hole ladder dia- 
grams by summing the combined particle-particle and hole- 
hole ladders to all orders in kpa while keeping the criteria (i) 
and (ii) satisfied. Such a resummation scheme for the unpo- 
larized case x — has been performed by Kaiser lf30tl - 

Following the treatment by Kaiser f50ll . we rewrite the 
propagator (fT9b in an alternative form 

GApo, P) = QoiPo, P) + 2mS(p - w p )0(A£ - |p|), (56) 

where the first term corresponds to the vacuum propagator 
£?o(po,p) = (po - w p + ie) and the second term is a so- 
called "medium-insertion" (MI) ll50ll . The elementary bubbles 
in this treatment are shown in Fig. [7] The first diagram Bo is 
identical to the vacuum part studied in Section III and it can 
be renormalized using the PDS scheme. For our purpose of 
resummation, we are interested in the following two quanti- 
ties, Bo + B\ + Z?2 and Bo + B\, which are mutually complex 
conjugate. We have 

Mix Mkp 

Bo + B x +B 2 = --r- + —±[R{ s j)-inI(s,fj\, 
4n 4n z 

Bo + Bi = -^tL + ^L [ R ( S , t) + frr/( Jt f)] t (57 ) 
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where R(s, t) = R-[(s, t) + Ry(s, t) and I(s, t) is the function 
defined in (30). 

To sum all ladder diagrams built from the elementary bub- 
bles, we first notice that the non-vanishing contributions to the 
interaction energy come from diagrams with at least two ad- 
jacent Mis l50h . Then a typical n-th-order contribution would 
look like the ring diagram of Fig. [TJ(b) with n vertices and at 
least two adjacent Mis. Naively, all these n-th-order diagrams 
are summed to give g"[(B + B\ + B 2 )" - (B + B{f] where 
the subtraction gets rid of those diagrams which have no adja- 
cent MI pairs. However, this expression is complex and there- 
fore cannot be the correct one. The crucial observations are 
that: (1) Each n-th-order ring diagram has a zi-rotational sym- 
metry. Therefore, we should introduce an additional factor 
1/n; (2) An n-th-order ring diagram comes from closing two 
open Mi-lines of an n-th-order ladder diagram, which intro- 
duces an integration over the allowed phase space |pi| < kl 
and |p2| < k F , but does not contribute a factor Bo to the energy 
as the naive expression does. These amendments lead to the 
correct n-th-order contribution to the interaction energy iBol : 
g n [(B () + Bi+ B 2 ) n - (Bo + Bi)"]/(2iln). The summation over 
n leads to two complex-conjugated logarithms and the final 
result is real. 

The final result for the energy density does not depend on 
the renormalization scale and the function f(x) in this re- 
summation scheme also takes the form (l43l , while the func- 
tion F(s, t) becomes 

In f 1 - -k F aR(s, t) + ik F aI(s, t)] - c.c. 

F(s, f) = -t n - — — . (58) 

2il(s, t) 

For small gas parameter k F a <k 1, F(s, t) can be expanded as 

1 , , 
F(s, t) = k F a + -(k F a) 2 R(s, t) + 0((k F aY). (59) 
n 



We can also check that 



so r , , r° 

—5- s as I 

n Jo Jo 



tdtl(s, t)R(s, t) = f^Ii^l , (60) 



which reflects the fact that the hole-hole ladders start to con- 
tribute at the order 0((k F a) 3 ) JM EH El • Therefore, the cri- 
teria (i) and (ii) are also fulfilled in the present resummation 
theory. Numerically, we also find a second order phase tran- 
sition, which occurs at a smaller gas parameter k F a = 0.786. 
We note that the inclusion of hole-hole ladders may not im- 
prove the quantitative result, since it only includes part of the 
beyond-leading-order contribution in the large-© expansion. 



VI. SUMMARY 

In summary, we have studied the non-perturbative effects 
on the ferromagnetic phase transition in repulsive Fermi gases 
by summing the ladder diagrams to all orders in the gas pa- 
rameter k F a. The non-perturbative effects not only reduce 
the critical gas parameter but also change the order of the 
phase transition. The resummation of particle-particle lad- 
ders, which corresponds to the leading order of the large- 
dimension expansion, predicts a second order phase transition 
occurring at k F a — 0.858, in good agreement with the Quan- 
tum Monte Carlo result 15J. The spin susceptibility calculated 
from our resummation theory are also in good agreement with 
the Quantum Monte Carlo results. Therefore, the resumma- 
tion of the ladder diagrams provides a more quantitative way 
to study the ferromagnetic transition in repulsive Fermi gases. 
In this paper, we have only considered a zero range poten- 
tial model. It will be interesting to study the non-universal 
shape-dependent contributions using the finite-density effec- 
tive range expansion ll47ll . 

Acknowledgments: We thank S. Pilati and S. -Y. Chang for 
providing us with the QMC data, N. Kaiser for helpful com- 
munications, and A. Sedrakian for reading the manuscript. 
LH acknowledges the support from the Alexander von Hum- 
boldt Foundation, and XGH is supported by the Deutsche 
Forschungsgemeinschaft (Grant SE 1836/1-2). 



[1] E. Stoner, Phil. Mag. 15, 1018 (1933). 

[2] K. Huang, Statistical Mechanics, Wiley, New York, 1987. 

[3] G.-B. Jo, Y.-R. Lee, J.-H. Choi, C. A. Christensen, T. H. Kim, 

J. H. Thywissen, D. E. Pritchard, and W. Ketterle, Science 325, 

1521 (2009). 

[4] C. Sanner, E. J. Su, W. Huang, A. Keshet, J. Gillen, and W. 

Ketterle, arXiv: 1108.2017. 
[5] S. Pilati, G. Bertaina, S. Giorgini, and M. Troyer, Phys. Rev. 

Lett. 105, 030405 (2010). 
[6] S.-Y. Chang, M. Randeria, and N. Trivedi, Proc. Natl. Acad. 

Sci. 108, 51 (2011). 
[7] L. Salasnich, B. Pozzi, A. Parola, and L. Reatto, J. Phys. B33, 

3943 (2000). 

[8] T. Sogo and H. Yabu, Phys. Rev A66, 04361 1(2002). 
[9] G. J. Conduit and B. D. Simons, Phys. Rev. Lett. 103, 
200403(2009). 

[10] J. L. LeBlanc, J. H. Thywissen, A. A. Burkov, and A. 



Paramekanti, Phys. Rev. A80, 013607 (2009). 
[11] S. Zhang, H.-H. Hung, and C. Wu, Phys. Rev. A82, 053618 
(2010). 

[12] H. Dong, H. Hu, X.-J. Liu, and P. D. Drummond, Phys. Rev. 

A82, 013627 (2010). 
[13] M. Okumura, S. Yamada, M. Machida, H. Aoki, Phys. Rev. 

A83, 031606(R) (2011). 
[14] R. A. Duine and A. H. MacDonald, Phys. Rev. Lett. 95, 230403 

(2005). 

[15] G. J. Conduit and B. D. Simons, Phys. Rev. A79, 053606 
(2009); G. J. Conduit, A. G. Green, and B. D. Simons, Phys. 
Rev. Lett. 103, 207201 (2009). 

[16] B. M. Fregoso and E. Fradkin, Phys. Rev. Lett. 103, 205301 

(2009) . 

[17] X.-J. Liu, H. Hu, and P. D. Drummond, Phys. Rev. A82, 023619 

(2010) . 

[18] H. Heiselberg, Phys. Rev. A83, 053635 (2011). 



12 



[19] H. Zhai, Phys. Rev. A80, 051605(R) (2009); X. Cui and H. 

Zhai, Phys. Rev. A81, 041602(R) (2010). 
[20] S. Q. Zhou, D. M. Ceperley, and S. Zhang, Phys. Rev. A84, 

013625 (2011). 

[21] D. Pekker, M. Babadi, R. Sensarma, N. Zinner, L. Pollet, M. W. 

Zwierlein, and E. Dernier, Phys. Rev. Lett. 106, 050402 (2011). 
[22] I. Sodemann, D. A. Pesin, and A. H. MacDonald, Phys. Rev. 

A85, 033628 (2012). 
[23] A. Recati and S. Stringari, Phys. Rev. Lett. 106, 080402 (2011). 
[24] G. J. Conduit and E. Altaian, Phys. Rev. A83, 043618 (201 1). 
[25] N. D. Drummond, N. R. Cooper, R. J. Needs, and G. V. Shlyap- 

nikov, Phys. Rev. B83, 195429 (2011). 
[26] X.-J. Liu and H. Hu, Phys. Rev. A82, 043626 (2010). 
[27] C. W. von Keyserlingk and G. J. Conduit, Phys. Rev. A83, 

053625 (2011). 

[28] C.-C. Chang, S. Zhang, and D. M. Ceperley, Phys. Rev. A82, 

061603(R) (2010). 
[29] G. Carleo, S. Moroni, F. Becca, and S. Baroni, Phys. Rev. B83, 

06041 1(R) (2011). 
[30] G. J. Conduit, Phys. Rev. A82, 043604 (2010). 
[31] R. A. Duine, M. Polini, H. T. C. Stoof, and G. Vignale, Phys. 

Rev. Lett. 104, 220403 (2010). 
[32] E. Taylor, S. Zhang, W. Schneider, and M. Randeria, Phys. Rev. 

A84, 063622 (2011). 
[33] F. Arias de Saavedra, F. Mazzanti, J. Boronat, and A. Polls, 

Phys. Rev. A85, 033615 (2012). 



[34] P. Massignan and G. M. Bruun, Eur. Phys. J. D65, 83 (201 1). 
[35] M. Sandri, A. Minguzzi, and F. Toigo, EPL 96, 66004 (2011). 
[36] F. Palestini, P. Pieri, and G. C. Strinati, Phys. Rev. Lett. 108, 
080401 (2012). 

[37] D. Belitz, T. R. Kirkpatrick, and T. Vojta, Phys. Rev. Lett. 82, 

4707 (1999). 

[38] J. V. Steele, arXiv: |nucl-th/00 10066] 

[39] S. Kanno, Prog. Theor. Phys. 44, 813 (1970). 

[40] K. Huang and C. N. Yang, Phys. Rev. 105, 767 (1957); T. D. 

Lee and C. N. Yang, Phys. Rev. 105, 1119 (1957). 
[41] H. W. Hammer and R. 1. Furnstahl, Nucl. Phys. A678, 277 

(2000). 

[42] D. B. Kaplan, M. J. Savage, and M. B. Wise, Nucl. Phys. B534, 
329 (1998); Phys. Lett. B424, 390 (1998). 

[43] E. Braaten and H. W. Hammer, Phys. Rept. 428, 259 (2006). 

[44] L. Pricoupenko and Y. Castin, Phys. Rev. A69, 051601 (2004). 

[45] A. L. Fetter and J. D. Walecka, Quantum Theory of Many- 
Particle Systems, McGraw-Hill, New York, 1971. 

[46] H. Heiselberg, Phys. Rev. A63, 043606 (2001). 

[47] T. Schafer, C.-W. Kao, and S. R. Cotanch, Nucl. Phys. A762, 
82 (2005). 

[48] V. B. Shenoy and T.-L. Ho, Phys. Rev. Lett. 107, 210401 

(2011). 

[49] P. Niemann and H.-W. Hammer, larXiv: 1 203 . 1 8241 
[50] N. Kaiser, Nucl. Phys. A860, 41 (2011). 



